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Abstract 

We present a review of atmospheric muon flux and energy spectrum measurements 
over almost six decades of muon momentum. Sea-level and underground/water/ice ex- 
periments are considered. Possible sources of systematic errors in the measurements 
are examinated. The characteristics of underground/water muons (muons in bundle, 
lateral distribution, energy spectrum) are discussed. The connection between the at- 
mospheric muon and neutrino measurements are also reported. 
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1 Introduction 

Primary cosmic rays (CRs) are particles accelerated at astrophysical sources continuously 
bombarding the Earth. When entering the Earth's atmosphere, primary CRs interact with 
the air nuclei and produce fluxes of secondary, tertiary (and so on) particles. All these par- 
ticles together create a cascade, called air shower. As the cascade develops longitudinally 
the particles are less and less energetic since the energy of the incoming CR is split and 
redistributed among more and more participants. The transverse momenta acquired by the 
secondaries cause the particles to spread laterally as they propagate in the target. Muons 
predominantly originate from the decay of secondary charged pions and kaons. The most 
important decay channels and respective branching ratios are: 

tt* -> fi ±( ul (~ 100%) (la) 

K± -»• n ±{ vl (~ 63.5%) (lb) 

Atmospheric muons are the most abundant charged particles arriving at sea level and 
the only ones able to penetrate deeply underground. The reason relies on their small energy 
loss in the whole atmosphere (~ 2 GeV), the relatively long lifetime and the fairly small 
interaction cross section. Because of the close relation between muon and neutrino pro- 
duction, the parameters characterizing muon physics can provide important information on 
atmospheric neutrino flux. 

An important parameter to describe the interactions and the subsequent propagation 
of the particles produced is the atmospheric depth X, measured in g/cm 2 , defined as the 
integral in altitude of the atmospheric density above the observation level h: 

POO 

X= / p(ti)dti ~ X e- h/ho (2) 

Jh 

In the last step, an approximation for an isothermal atmosphere was used, where Xo = 
1030 g/cm 2 is the atmospheric depth at sea level and ho — 8.4 km is the scale height 
(for mid latitudes) in the atmosphere. Eq. ([2]) is valid for vertically incident particles. 
For zenith angles < 60°, for which the Earth surface can be approximated as flat, the 
atmospheric depth is scaled with 1 / cos 8, giving the slant depth. For larger zenith angles, 
the curvature of the Earth has to be accounted for. The atmospheric profile gives a total 
horizontal atmospheric depth of about 36000 g/cm 2 . 

The air shower is described by a set of coupled cascade equations with boundary condi- 
tions at the top of the atmosphere to match the primary spectrum. Using the transport equa- 
tions, analytic expressions of the cascade can be constructed. The solutions of these equa- 
tions allow to compute the differential particle flux anywhere within the atmospheric target. 
Some approximate analytic solutions are valid in the limit of high energies [39], [40], [50]. 
Numerical or Monte Carlo calculations are needed to account accurately for decay and en- 
ergy loss processes, and for the energy-dependences of the cross sections and of the primary 
spectral index. The nucleon mean free path Ajv in atmosphere is given (in g/cm 2 units) by 

Am n 

= — - (3) 

a N 

where a^ r is the interaction cross section of nucleon in air, A is the mean mass number 
of air nuclei and m p the proton mass. For nucleons in the TeV range, a^ r ~ 300 mb. In 
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the context of air shower development, the energy-dependent cross section for an inelastic 
collision of a nucleon with an air nucleus is assumed to be constant. The atmosphere of 
the Earth consists mainly of nitrogen and oxygen: the interaction target for the primary 
beam is half protons and half neutrons. Assuming an average atmospheric nucleus with 
A ~ 14.5, Atv — 80 g/cm 2 . The total vertical atmospheric depth is about 1000 g/cm 2 and 
it corresponds to more than 1 1 interaction lengths. 

Most muons are produced through processes ([TJ high in the atmosphere in the first few 
interaction lengths. The decay mean free path of pions d n in units of slant depth, is defined 
as: 

J_ _ m n c 2 h _ e n ^ 
d n Ecr n X cos 9 EX cos 6 

where E, , t w are the pion energy, mass and lifetime, respectively. A similar relation 
holds for the kaon. Decay or interaction dominates depending on whether l/d w or 1/X n 
is larger. X w is defined through Eq. ([3]) replacing <tjv with a n . At the critical energy 
E = e n = m n c 2 ho/T 7T c = 115 GeV the interaction probability in the atmosphere equals 
the decay probability. 

As for the pion, all long-lived unstable particles [24] are subject to competition between 
interaction and decay as they propagate in the atmosphere. The probability for either pro- 
cess to occur depends on the lifetime of the particle and is a function of its kinetic energy 
and on the local atmospheric density, which is a function of altitude. This interrelationship 
is responsible for the zenith angle enhancement of the bulk of the muons in air showers. 

The zenith angle enhancement phenomenon does not affect the distribution of muons 
produced in semileptonic decays of charmed mesons, like D^, D° and others. As the life- 
time of charmed particles is smaller than ~ 10 -12 s {prompt decays), they yield so-called 
prompt (or direct) muons that are in general highly energetic for kinematic reasons. Since 
the production cross section of charmed mesons in proton-nucleon interactions is rather 
small, D decays contribute significantly only at very high energies. 



2 Cosmic rays at the sea level 

Muons are the dominant component of charged particles at sea level. The integral fluxes 
of particles arriving at geomagnetic latitudes ~ 40° vs. their kinetic energy are presented 
in Fig. [T] Fluxes are averaged over the 11 -year solar cycles. The muon flux with E^ > 1 
GeV through a horizontal area amounts to roughly one particle per cm 2 and per minute: 
I v {En > 1 GeV) ~ 70 m~ 2 s" 1 sr" 1 E21 

The muon energy spectrum at sea level is a direct consequence of the meson source 
spectrum. Below the few GeV energy range, the muon decay probability cannot be ne- 
glected. A muon of 1 GeV has a Lorentz factor 7 = E^/m^c ~ 10 and it has a mean 
decay length d^ = jt^c ~ 6 km. Since pions are typically produced at altitudes of 15 
km and decay relatively fast (for 7 = 10 the decay length is d w ~ 78 m, which is al- 
most the same value of X w ), the daughter muons do not reach the sea level but rather decay 
themselves or get absorbed in the atmosphere. 

At higher energies, the situation changes. For pions of 100 GeV (d n ~ 5.6 km, cor- 
responding to a column density of 160 g/cm 2 measured from the production altitude) the 
interaction probability starts to dominate over decay. Pions of these energies will therefore 
produce further tertiary pions in subsequent interactions, which will also decay eventually 
into muons, but providing muons of lower energy. Therefore, the muon spectrum at high 
energies is always steeper compared to the parent pion spectrum. 
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Figure 1: Integral fluxes averaged over the 11 -year solar cycle of fj,, e,p and photons (ph) 
arriving at geomagnetic latitudes ~ 40° vs. their kinetic energy. 

The muon energy and angular distribution is the effect of a convolution of production 
spectrum, energy losses in the atmosphere and decay. The competition of decay and in- 
teraction plays a crucial role and the relative importance of the two processes depends on 
energy. The mean energy of muons at the ground is about 4 GeV 11511 . 

Three different energy regions in the muon spectrum are distinguishable: 

• < e M , where e M ~ 1 GeV. In this case, muon decay and muon energy loss are 
important and must be taken into account. Semi-analytical solutions (as those used in the 
higher energy range) overestimate the flux. The energy spectrum is almost flat, starting to 
steepen gradually in the same way as the primary spectrum above 10 GeV. 

• C/i < Efj, < e-w^K, where e^= 115 GeV and ex = 850 GeV are the critical energies for 
the vertical directions. In this energy range, almost all mesons decay and the muon flux has 
the same power law of the parent mesons, and hence of the primary CRs. The muon flux is 
almost independent on the zenith angle. 

• En 3> e n : k. The meson production spectrum has the same power law dependence of the 
primary CRs, but the rate of their decay steepen one power of E^ since the pion and kaon 
decay probability is suppressed. The thickness of the atmosphere is not large enough for 
pions to decay, since the high Lorentz factor. For E > the inclined muon spectrum is 
flatter than the vertical one and the muon flux is respectively higher. 

In the intermediate and high energy region (above 100 GeV) and for zenith angle 6 < 
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60° an approximate formula holds: 



dE a dQ 




l + (^)cos0 l + (^)cosfl 



1 B 



(5) 



where the scale factor A, the power index 7, the balance factor B (which depends on the 
ratio of muons produced by kaons and pions), and the a, b coefficients are adjustable param- 
eters. Different best estimates of these parameters were published by several authors. For a 
review, see ||49| . Because pions decay more easily in non- vertical showers, a zenith angle 9 
factor enters in the formula and muons at large angles have a flatter energy spectrum. 

2.1 Angular distribution at sea level 

The muon intensity from horizontal directions at low energies is naturally reduced because 
of muon decays and absorption effects in the thicker atmosphere at large zenith angles. At 
high energy the parent particles of muons travel relatively long distances in rare parts of 
the atmosphere. As a consequence, their decay probability is increased compared to the 
interaction probability. 

Fig. [2] gives a quantitative description of this effect. Muons below the few GeV/c mo- 
mentum range fade fairly quickly with increasing zenith angle, with dependence oc cos n 9, 
where n ~ 2 3. 

The flux of muons in the 100 GeV/c range is relatively flat up to cos# ~ 0.2 and 
then quickly declines. At 1 TeV/c the flux monotonically increases with the zenith angle, 
approaching the 1 / cos 9 dependence. The flux of TeV muons is particularly sensitive to 
large values of the zenith angle. When approaching the horizontal direction, a small differ- 
ence in cos 9 changes appreciably the thickness and the density profile of the atmosphere 
and the corresponding muon energy spectrum. For this reason the measurements of almost 
horizontal muons is very difficult. 

3 Measurements at sea level 

Measurements performed at ground level offer the advantage of a high stability, large col- 
lecting factor and long exposure time due to relatively favourable experimental conditions. 
Sea level data offer the possibility to perform a robust check of the reliability of existing 
Monte-Carlo codes. 

Many of the experiments devoted to the measurement of the muon momentum spectra 
and intensity (vertical and inclined directions) have been carried out since the '70s. The 
results are often in disagreement with one another; the discrepancies are significantly larger 
that the experimental reported errors. Recently new instruments, mainly spectrometers de- 
signed for balloon experiments or used primarily in CERN LEP and LHC experiments and 
used also to CR studies, have added new valuable information. 

The vertical muon intensity at sea level is a quantity that varies with the geomagnetic 
latitude, altitude, solar activity and atmospheric conditions. When comparing muon obser- 
vations at low energies (< 20 GeV/c) it is important to know the year and location where 
the measurements were made ll27ll . 

The geomagnetic field tends to prevent low energy cosmic rays from penetrating 
through the magnetosphere down to the Earth's atmosphere. At any point on the Earth 
one can define a threshold (or cut-off) rigidity for cosmic rays arriving at a particular zenith 
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Figure 2: Angular distribution of muons at the ground for different muon energies. The 
overall angular distribution of muons measured at sea level is oc cos 2 9, which is charac- 
teristic of muons with Ep ~ 3 GeV. At lower energy the angular distribution becomes 
increasingly steep, while at higher energy it flattens, approaching a 1/ cos 9 distribution for 
En S> e-w and 9 < 60°. At large angles low energy muons decay before reaching the surface 
and high energy pions decay before they interact, thus the average muon energy increases. 
The estimate of the angular distribution is based on a Monte Carlo and accounts for the 
curvature of the Earth atmosphere. 

and azimuth angle. Primary nuclei having lower rigidity are deflected by the action of the 
geomagnetic field and do not produce muons at that latitude. The cut-off values range from 
less than 1 GV near the geomagnetic poles to about 16 GV for vertical particles near the 
equator. As CR primaries are predominantly protons and nuclei, it results that at a given 
location the intensity from the West is stronger than that from the East. The geomagnetic 
effects are important for sea level muons up to about E^ ~ 5 GeV, and the effect is larger 
at higher altitudes. 

The 1 1 year solar cycle influences the primary CR spectrum at the top of the atmo- 
sphere, as the configuration of the Interplanetary Magnetic Field varies. It results that the 
cosmic ray flux is significantly modulated up to energies of about 20 GeV. 

Most experiments are not exactly performed at sea level. A correction to take into ac- 
count the dependence of the flux from the altitude must be included. For muon momenta 
above 10 GeV and altitudes H less than ~ 1000 m the vertical muon flux can be parame- 
terized as I^H) = Ip(0)e- H / L ( p \ where L(p) = 4900 + 750p is a scale factor (in meters) 
which depends on muon momentum, p (measured in GeV/c). 



Atmospheric muons: experimental aspects 



7 



1E-02 



" 1E " 03 



1E-04 



1E-05 



30 



A Flint 
XJakeman 
■ Wilson 

♦ Gettert 

* Dmitrieva 

• Cro&Ras 



40 



50 60 70 

zenith angle [°] 



SO 



Figure 3: Variation of absolute integral intensity of muons at sea level with the zenith 
angle. The momentum threshold pr is slightly different for each experiment and ranges 
from 0.1 to 1 GeV/c. The points are normalized assuming px = 0.35 GeV/c. The points 
correspond (from top to bottom): l38l. fi7l 051. fiTI 051 



Changes in pressure and temperature in the atmosphere above the detector produce 



variations which will be considered in i 5.2 



3.1 Experimental setups 

Different experimental methods were used to measure the muon flux and energy spectrum. 

Muon telescopes are made of several charged particle detectors arranged along a straight 
line and interlaid by one or more layers of absorbing material. In some experiments the de- 
tector and absorber are in a rigid construction which could be rotated in zenith and azimuth 
direction, allowing the selection of muons from a given direction of celestial hemisphere. 
The quantity of material (in g/cm 2 ) travelled by muons in such a telescope is approximately 
constant and it sets the muon energy threshold. 

Multi-directional muon telescopes generally consist of at least two layers of segmented 
muon detectors. The coincidence of signals between two counters in upper and bottom 
layers determines the direction of muon arrival. The quantity of material crossed by the 
particle in such detectors increases with increase of zenith angle, so the threshold energy 
for multidirectional muon telescopes depends on 9. A compilation of measurements of the 
muon flux at latitudes between 52° and 56° as a function of the zenith angle by different 
experiments is reported in Fig. [3] 

The muon energy spectrum has been extensively measured, mainly by solid iron mag- 
net spectrometers. For these detectors, the multiple scattering plays an important role in the 
momentum resolution, particularly at low energies. Recently, measurements have been per- 
formed using low mass superconducting magnet spectrometers designed as a balloon-borne 
apparatus for cosmic ray studies. They represent a sort of second generation experiments. 

The atmospheric muon flux and energy spectrum was also measured using the precise 
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Reference 


Vll [GeV/c] 


B 


H&T 


PDG 


[13] 


15.1 - 82.1 


X 


0.79 




El 


20-500 


X 




X 


El 


11-810 


X 


X 




ii 


3- 3000 


X 


0.933 


X 


El 


10-150 


X 


0.858 







20-1000 


X 


1.039 


X 


[32]* 


0.25-100 


X 


0.944 


X 




0.2-120 




0.818 


X 


h 


20-3000 






X 


B5I* 


0.6-400 






X 



Table 1 : Compilation of different measurements of the muon momentum spectrum. Ex- 
periments with magnetic spectrometers are indicated with *. The second column reports 
the momentum range; the 3rd,4th and 5th columns if the data are used in the compilation 
of 11231 , P31l (with the used normalizations factor) and lIBTI , respectively. 



muon spectrometer of the L3 detector which was located at the LEP collider at CERN. This 
apparatus collected muons ~30 m below a stratified rock overburden, and with dimensions 
much larger than previous experiments (the volume of the region with a magnetic field of 
0.5 T was ~ 1000 m 3 ). 

Finally, information about the muon flux with momenta larger than few TeV have been 
extracted from underground measurements, see £j4j Comprehensive review of various types 
of particle detectors used for cosmic ray studies can be found in ||37l . 

3.2 Momentum spectra at the vertical direction 

Table [T] lists the authors and energy range of the reported absolute vertical intensity mea- 
surements. It is also shown if the experiment has been used by other reviews, namely "B" 
11231 . "H&T" [45] and in the Particle Data Group "PDG" HQ. For H&T we report also 
the final normalization factor they have found. The symbol (*) refers to experiment using 
superconducting magnet spectrometers. 

The measurements listed in Table Q] (only published results) of the muon momentum 
from the vertical direction are presentea in Fig. [4] The agreement between measurements 
is relatively good and the largest contribution to the deviations are the systematic errors 
due to incorrect knowledge of the acceptance, efficiency of the counters and corrections for 
multiple scattering. 

Measurements of the muon momentum spectra for p M < 1 TeV/c are particularly im- 
portant for the comparison of nuclear cascade models with available data. The differences 
between the measurements of the sea-level spectra are more evident if the percentage de- 
viations of the data from one of the parameterization are plotted (see Fig. [5]). Differences 
of ± 15-20% exist between the data and the parameterization. The disagreement between 
the different experiments can be as large as 30-35% despite the fact that the declared indi- 
vidual errors are small (however increasing with momentum due to the decreasing number 
of detectable particles and to the maximum detectable momentum). The origin of the dis- 
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Figure 4: Vertical absolute differential muon intensity p^dN^/dE^dQ at sea level for the 
experiments reported in Tab. [I] The ordinate values have been multiplied by p^ in order to 
compress the plot and to emphasize the differences. In this energy range, p M ~ E^. 



crepancies can well be due to the incorrect knowledge/control of systematic effects. For 
example, in the case of L3+C the stable negative deviations can be due to a bad correction 
for the molasses that cover the apparatus. 



4 Underground muons 



Underground measurements offer the possibility to extend the energy range of the muon 
spectrum beyond 1 TeV. Such data are of an indirect type, but their link with the direct 
lower energy observations gives the possibility to complete the picture of muon spectra 
measurements. 

Deep underground detectors have normally large collecting area and are not subject to 
the time restrictions of balloon and satellite experiments, so they can measure the muon 
flux for a very long time. As a drawback, as discussed in §2, the energy dependence of 
the muon spectrum at E^ > 1 TeV steepen one power and the intensity is a factor 1000 
smaller than that of primary CRs on the upper atmosphere. The penetrating component of 
CRs underground depends on a complex convolution of different physics processes, as the 
high energy muon production spectrum and the muon energy losses. Particularly important 
is the knowledge of the composition and thickness of the material overburden the detector. 
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Figure 5: Relative deviations of the differential muon intensity measurements at sea level 
with respect to the [ 23 ] parameterization. The symbols are the same as in Fig. [4j the indirect 
measurements described in £4.2 and Fig. [7] are also included. 



4.1 Muon energy losses 

Muon energy losses are usually divided into continuous and discrete processes. The former 
is due to ionization, which depends weakly on muon energy and can be considered nearly 
constant for relativistic particles. For muons below ~ 500 GeV, this is the dominant energy 
loss process. For muons reaching great depths, discrete energy losses become important: 
bremsstrahlung (br), direct electron-positron pair production (pair) and electromagnetic 
interaction with nuclei (photoproduction, ph). In these radiative processes energy is lost in 
bursts along the muon path. In general the total muon energy loss is parameterized as: 

(6) 

where X is the thickness of crossed material in g/cm 2 and /3 = Pbr+fipair+Pph is the sum of 
fractional energy loss in the three mentioned radiation processes. As the rock compositions 
are different for different underground experiments, the so-called standard rock is defined 
as a common reference. The standard rock is characterized by density p = 2.65 g/cm 3 , 
atomic mass A = 22 and charge Z = 1 1 . The thickness X is commonly given in units of 
meters of water equivalent (1 m.w.e.= 10 2 g/cm 2 ). 

The factors a and j3 in Eq. ([6]) are mildly energy dependent as well as dependent upon 
the chemical composition of the medium: in particular a oc Z/A and oc Z 2 /A. Typical 
values are a ~ 2 MeV g _1 cm 2 and /3 ~ 4 x 10~ 6 g _1 cm 2 . The critical energy is defined as 
the energy at which ionization energy loss equals radiative energy losses: e M = a/f3 ~ 500 
GeV. 
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Figure 6: Minimum energy at surface, E m i n to reach a given underground depth X = h. 
Typical values of E m i n for muons are 2, 20 and 2000 GeV to reach the ground (10 m.w.e.), 
the LEP tunnel (~80 m.w.e.) and the Gran Sasso Laboratory (~3 km.w.e.), respectively. 



The general solution of Eq. (|6j) corresponds to the average energy (E^) of a beam of 
muons with initial energy E® after penetrating a depth X: 



(E lx (X)) = (E°+e,)e 



-I3X 



(7) 



The minimum energy required for a muon at the surface to reach slant depth X is the 
solution of Eq. ^ with residual energy E^(X) = 0: 



E 



ex 



(8) 



The range R for a muon of energy E®, i.e. the underground depth that this muon will reach, 



is: 



1 K 



(9) 



The above quantities are average values. For precise calculations of the underground muon 
flux one needs to take into account the fluctuations inherent to the radiative processes. Be- 
cause of the stochastic character of muon interaction processes with large energy transfers 
(e.g., bremsstrahlung) muons are subject to a considerable range straggling. The higher E® 
is, the more dominant are the radiation processes and the more important are the fluctua- 
tions of the energy losses which broaden the distribution of the range. Fig. [6] shows the 
typical values of the minimum energy at surface, E m i n to reach a given underground depth. 
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Experiment 


Depth 


Momentum 


(Reference) 


(m.w.e.) 


(TeV/c) 


LVD Q 


>3000 


1.9-43 


MACRO QQ 


>3150 


0.5-20 


Baksan OH 


>850 


1-30 


MSU 1551 


>50 


3-50 



Table 2: Compilation of different indirect measurements of the muon momentum with 
underground detectors. 



4.2 The deep-intensity relation 

The muon spectrum at energies above few TeV was deduced by underground measure- 
ments. The procedure used for this indirect measurement of the sea-level energy spectrum 
passes through the determination of the so-called deep-intensity relation (DIR) I®(h). This 
quantity represents the muon intensity at the vertical direction as a function of the depth h. 
As underground detectors are at a fixed position, in principle only one point can be mea- 
sured. However, when measuring the muon intensity I^h, 9) at different zenith angle 9, 
the quantity of rock (or water) overburden changes. At high energy (E^ > 1 TeV) and for 
9 < 60°, Eq. (|5j> provides a simple relationship between I®(h) and I^h, 9): 

I (t (h,9)=I°(h)/coB0 (10) 



From the experimental point of view, I^(h, 9) is measured at a given direction 6, corre- 
sponding to a slant depth h, as: 

t ( h n\ = ( 1 \ Ei N i m i nn 

V AT / Y.j AfJj(0, 4>)Aj(9, (j))ej(9, (ft) { 1 

where AT in the total livetime (in s) of the experiment, Ni is the number of detected events 
with multiplicity mi with zenith 9 and azimuth <p in the angular bin Aflj (9, <ft) (sr). Aj (9, (ft) 
(m 2 ) and ej(9, (ft) correspond to the geometrical intrinsic acceptance of the detector, and the 
overall efficiency in the (9, (ft) bin. 

The relation between the measured DIR function I^{h) (measured from Eq. (Ill with 

the use Eq. §) and the differential sea-level muon spectrum ™fc is expressed^ 

f 00 dN 

'> )= L i^n p(E " h)dE '- (12) 



Here, P(E^, h) is the muon survival function. It represents the probability that muons of 
energy at surface reach a given depth h, and it is determined via Monte Carlo calcula- 
tions. Assuming Eq. (|5]> for the sea-level muon spectrum, leaving as free parameters the 
muon spectral index 7, the scale A and the balance B factors, it is possible to unfold the 
sea level muon spectrum from the measured vertical muon intensity. 
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Figure 7: Extension of the data shown in Fig. [4] to higher momenta, by including the 
indirect measurement of the underground muon flux. 

Table [2] indicates the experiments, the depth and the estimated muon momentum range 
for these indirect underground measurements. The results are shown in Fig. [7] The differ- 
ential distribution (dN^/dE^dQ) is presented, as usual, multiplied by a factor p 3 ^ (as mo- 
mentum and energy coincide) to better observe the variation of the spectrum in the whole 
energy region. In these indirect measurements, the main sources of systematic uncertainties 
are due to the treatment of hard processes in the energy loss of muons and to the knowledge 
of the rock density and overburden thickness, which rely on geological surveys. 

This situation seems to favour underwater/ice experiments. In this case, the uncertainty 
in the material density and overburden is negligible with respect to underground experi- 
ments. However, when small-size experiments are carried out, the main uncertainty is the 
multiplicity of the muon bundle. In the case of large underwater experiments (like the 
running neutrino telescopes, (281), they are optimized to look for upward-going neutrino 
induced particles. Atmospheric muons are seen with the "tail of the eyes" (the photomul- 
tipliers inside optical modules), where large uncertainties on the optical module angular 
acceptance do not allow a precision measurements. See [6] for a discussion and the DIR 
measurement with the ANTARES underwater neutrino telescope. 



5 Characteristic of underground/underwater muons 



Fig. [8] shows a multiple muon event detected by MACRO. Multiple events are closely 
packea bundles of parallel muons, usually of high energy, originating from the same pri- 
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Figure 8: A multimuon event seen in the MACRO experiment. 10 different tracks are 
identified. 

mary CR. Multi-muon events are used to explore the properties of very high energy hadronic 
interactions and to study the longitudinal development of showers 0. The multiplicity of 
produced secondary particles increases with the energy of the initiating particle. The muon 
multiplicity is an observable which is also correlated with the mass of the primary CR: at a 
given total energy, heavier nuclei produce more muons than a primary proton. 

The interaction vertex of the particles which initiate the air showers is typically at an 
atmospheric altitude of 15 km. Since secondary particles in hadronic cascades have small 
transverse momenta p t (~ 300 MeV/c), high energy muons are essentially collimated near 
the shower axis. Considering a primary nucleon, producing mesons of energy E^jc with 
transverse momentum pt at a height H pro d, the average separation of their daughter high 
energy muons from the shower axis is given by ifTTTl : 

r ^ -^—H-prod (13) 

For primary energies around 10 14 eV, the lateral displacements of energetic muons (~ 1 
TeV) of several meters are typically obtained underground. Displacements are almost ex- 
clusively caused by transferred transverse momentum in hadronic processes. Typical mul- 
tiple scattering angles for muon energies around 100 GeV in thick layers of rock (50-100 
m) are on the order of a few mrad. 



5.1 Atmospheric muons in neutrino telescopes 

Atmospheric muons represent the most abundant signal in a neutrino telescope and can be 
used to calibrate the detector and to check its expected response to the passage of charged 
particles. On the other side, they can represent a dangerous background source because 
downward-going muons can incorrectly be reconstructed as upward-going particles and 
mimic high energy neutrino interactions; muons in bundles are particularly dangerous. 
These muons are expected to arrive almost at the same time in the plane perpendicular 
to the shower axis. A full Monte Carlo simulation, starting from the simulation of atmo- 
spheric showers, can accurately reproduce the main features of muons reaching a neutrino 
telescope, but requires a large amount of CPU time. 

Recently, parametric formulae lf20l to evaluate the flux of atmospheric muons were de- 
rived from a full Monte Carlo simulation. These formulas take into account the muon mul- 
tiplicity and the energy spectrum of muons in a bundle, as a function of the distance from 
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Figure 9: Flux of bundles of muons as a function of their multiplicity m in the bundle 



obtained with Eq. ( 14 1 at the depth of 3 km.w.e. and for 5 different values of the zenith 



angle. The computation assumes here a muon energy threshold of 1 GeV. 

the shower axis. A simple generator interface is provided (called MUPAGE), which can 
be used by all experiments having a flat overburden coverage of at least 1500 m.w.e. Eoll . 
MUPAGE is used to simulate atmospheric muons in Mediterranean neutrino telescopes. 

The flux (which corresponds to the (dN^/dE^dil) of Eq. ^ integrated over the muon 
energy E^, with units: m~ 2 s~ 1 sr~ 1 ) of muon bundles with multiplicity m (see Fig. 
is obtained as a function of the depth along the vertical direction ho (note that the index (J 
means that the depth h is computed exactly at the vertical direction) and zenith angle 6 as: 

m u(h, ,V) 

The flux of bundles of increasing multiplicity m decreases with increasing vertical depth 
and zenith angle. The parametric dependences of K(ho,0) and u(ho,0) are reported in 
ll20l . as the others below. 

The energy spectrum of muons is described IT391 by: 



5 ^- G ^™ [B( , + e( l-.-**) r (15, 
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Figure 10: Muon flux as a function of the cos 9 as measured by ANTARES0, AMANDA- 
II El, IceCube El, and Baikal EH at four different depths. The results of the MUPAGE 
computation is superimposed as full lines. 



where X = ho/ cos 6, 7 is the spectral index of the primary CRs and e = a//3; a, f3 are 
defined in Eq. ([6]). In the parameterization, 7 and e were instead considered as free fit 
parameters. The constant G = G(pj, e) represents a normalization factor, in order that the 
integral over the muon energy spectrum ( [T5] > from 1 GeV to 500 TeV is equal to 1. In the 
case on single muon events (bundles withm = 1) a simple dependence 7 = -y(ho) and 
e = e(/io, 9) holds. 

The situation is more complicated for multiple muons. Due to the muon production 
kinematics, the muon energy depends from their distance with respect to the axis of the 
bundle. The description of the muon lateral distance R from the axis is thus the preliminary 
step to evaluate the muon energy distribution in a bundle. R (in the plane orthogonal to the 
axis) was extracted from a distribution of the form: 



dN 
~d~R 



C 



R 



(r+Ro) c 



(16) 



where Rq = Ro(ho, m, 9) and a = a(/io, m). The energy spectrum of muons arriving in 
bundles has the same general form as for single muons ( 15 1. In the case of multiple muons, 
the analytic description of the parameters 7 = 7(/i, R,rn) and e = e(h, R, 9) depends on 
15 constants. 

These parameterizations allow to evaluate not only the total muon flux, but also the 
total number of muon bundles in deep detectors starting from the primary CR flux, CR 
composition and interaction model which reproduces (at the level of ~ 30%) the MACRO 
data (depth: 3000 -6000 m.w.e., 9 < 60°). Fig . [T0| shows the comparisons of the zenith 
distribution evaluated at a fixed depth using Eq. <\ 14|) and some underwater/ice data. 
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5.2 Pressure and temperature effect 

If atmospheric temperature changes by AT, the muon flux at the observation level X will 
change by a quantity AI M . This variation can be either positive or negative, i.e. AI^ oc 
±AT, depending on the energy of the observed muon. Two competitive effects are in 
action if the atmospheric temperature increases. As a consequence, the atmosphere expands 
and the air density decreases; the probability of the interaction of mesons (kaons and pions) 
at unit of geometric path becomes smaller, hence decay probability into muons becomes 
higher. On the other hand, the geometric expansion of the atmosphere increases the path 
from the generation point to the detector, and a higher number of muons will decay. For 
low energy muons, the latter effect is the dominant one, and the correlation sign between 
flux and temperature is negative (A/ M oc — AT). For high energies, muons have not enough 
time to decay in the atmosphere and the correlation sign becomes positive [36). 

Many under-ground/water/ice experiments measured the correlation between tempera- 
ture and high energy muon intensity. It was found that iTTOTl an effective temperature T e ff, 
defined by the weighted average of temperatures from the surface to the top of the atmo- 
sphere, is useful to describe the situation. T e f j approximates the atmosphere as an isother- 
mal body, weighting each pressure layer according to its relevance to muon production in 
atmosphere. The variation of muon rate AI^/I^ is related to the effective temperature as: 

= flT^ (17) 

where ay is the atmospheric temperature coefficient, which is a function of both the muon 
threshold energy and the K/tt ratio. As the energy increases, the muon intensity becomes 
more dependent on the meson critical energy e W: K, which in turn depends on the atmo- 
spheric temperature. The ar coefficient reflects the fraction of mesons that are sensitive to 
atmospheric temperature variations. For energies much greater than the critical energy, and 
thus for very deep experiments, the value of ax approaches unity. The expected effective 



temperature coefficient as a function of depth is shown in Fig. [TT] together with the values 
measured by underground/ice experiments. 



6 Muons and neutrinos 

As indicated in Eq. ([T]), the production mechanisms of atmospheric muons and neutrinos are 
strongly correlated, see [46] for a recent review. However, due to the two-body kinematics, 
the energy spectra of the /x's and z^'s from mesons decay are different. Let us consider 
for instance the pion decay in the centre of mass (cm.) system (m^ = 139.6 GeV; = 
105.7 GeV). The cm. muon energy is T* = (ml + mf 1 )/2m n = 109.8 MeV. Similarly 
for the neutrino, considering that in the cm. system E* + E* = m n , one has: E* = 
{m^ — m'j t )/2m 7T = 29.8 MeV. In the laboratory system, the energies are boosted by 
the Lorentz factor 7 = E^/m^. In any case, muons carry a larger fraction of the meson 
energy than neutrinos. As consequence, the energy distribution of is slightly shifted 
towards lower energy values than charged muons, as shown in Fig. fl2| Additional are 



produced by the in flight decay of muons, together with a u e and an electron/positron. As the 
muon decay probability in the atmosphere decreases with increasing E^, the u e spectrum is 
depleted with respect to that of at high energy. In Fig. [T2| we include the measurement 



of the f M energy spectrum reported by the Frejus, AMANDA and IceCube experiments 



18 



S. Cecchini and M. Spurio 



a T i 





























1 


1 

\ 


A ' 


*— 


¥■ 

* 






♦ MACRO 


• IceCube 


s 

s 

✓ i 


i 




m AMANDA XBaksan 
k MINOS FD • UTAH 
A LVD ASherman 

■ Borexino ■ Barrett 

l 


t 

/ 

< 


> 













1000 2000 3000 4000 

depth [m.w.e.] 



Figure 1 1 : The temperature coefficient ay as a function of detector depth. The dashed 
curve is the prediction using the pion-only model iTTOll . The points correspond (from top to 
bottom): iTTOl.HBI. E3, El, |(Bellini, 2012)1 , 01, H3, 61, El, EQ. 



The v e component was measured between (E^) ~0.4 - 14 GeV by the Frejus experiment. 
The measurement of the muon and neutrino energy spectra represents a very challenging 
result, as completely different experimental techniques were used to measure the charged 
and neutral leptons in different energy ranges. 

Neutrino telescopes are taking data under the Antarctic ice or in the Mediterranean sea. 
The main goals are the detection of neutrinos from cosmic sources and the measurement of 
the isotropic flux of high energy neutrinos from the ensemble of all extragalactic sources. 
The signature for the former is an excess of events over the background of atmospheric 
neutrinos from a given direction. The signatures of the diffuse astrophysical neutrino signal 
are: (a) isotropy; (b) a hard energy spectrum; (c) approximately equal fluxes of u e , and 
v T . The neutrino fluxes generated by the prompt charm decay have also the properties (a) 
and (b), and equal fluxes for v e and and therefore constitute a dangerous background. 



7 Conclusions 

The energy spectrum of atmospheric muons is measured over almost 6 decades of muon mo- 
mentum with different experimental techniques at sea level and in underground/water/ice 
experiments. The so-called prompt component, which is expected to dominate the spectrum 
at very high energies, is still undetected. Below 1 TeV/c, disagreements between different 
experiments are up to 20% due to systematic uncertainties. Slightly larger uncertainties 
arise from the indirect methods used to deconvolve the higher energy spectrum up to ~ 40 
TeV. These measurements of the energy spectrum above few TeV were made with large 
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underground detectors. In particular, the MACRO experiment at the Gran Sasso Laboratory 
in Italy accurately studied the multiplicity distribution of muons arriving in bundle, their 
lateral distribution with respect to the shower axis and their energy spectrum. These infor- 
mation were used to obtain a parameterization of the atmospheric muons detected by large 
underwater neutrino telescopes. 

The energy spectrum of atmospheric muon neutrinos are measured through the detec- 
tion of upgoing muons, generated by charged current interactions of atmospheric v^. They 
represent the irreducible background for searches of cosmic neutrinos, and for this reason 
this component must be accurately known. The u e component has still to be accurately 
measured, particularly in the high energy region. 

The knowledge of the atmospheric muon spectrum, the characteristics of the muon flux 
at sea level and the processes of attenuation of muons passing through rocks or other mate- 
rials is of fundamental importance for Earth science purposes and muon imaging feasibility. 
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Figure 12: Compilation of measurements of atmospheric muons and neutrinos. The points 
correspond (from top to bottom): E3, 0, Q, Q, ED, ED- The lines represent the 
neutrino flux expectation (full: vertical direction; dashed: horizontal direction) computed 

by era. 



